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ABSTRACT 

1.. 


A  class  of  models  for  target  patterns  (concentric  circular  waves 
emanating  from  a  point  called  the  leading  center)  is  constructed  in  the 
context  of  singularly  perturbed  reaction-diffusion  systems  of  partial 
differential  equations.  First,  the  theory  of  wave  fronts  is  detailed  for 
scalar  equations  and  systems  of  equations.  A  scaling  method  reduces  complex 
waves  to  the  consideration  of  a  group  of  simple  wave  phenomena.  It  is  shown 
that  expanding  wave  fronts  can  be  generated  spontaneously  at  a  point.  This 
process,  together  with  the  laws  of  their  subsequent  motion,  reduces  the 
problem  to  an  ordinary  differential  initial  value  problem,  whose  solution  is 
required  to  have  certain  properties.  A  discussion  is  given  of  the  connection 
between  these  results  and  experimental  observations. 
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•This  is  a  revised  and  expanded  version  of  a  paper  presented  by  Fife  at  a 
conference  with  proceedinqs  published  as  Dynamics  of  Synergetic  Systems,  H. 
Haken,  ed,  Springer-Verl  (1980). 
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SIGNIFICANCE  AND  EXPLANATION 


Propagation  chemical  waves,  particularly  those  enanating  from  a  "leadina 
center"  or  those  forminq  a  spiral,  may  he  seen  in  laboratory  reagents 

V*  v 

involving  the  Pelousov-Xhabotinskii  reaction.  Chemical  and  physico-chemical 
waves  also  occur  in  bioloaical  media.  Typically,  these  phenomena  exhibit 
multiple  natural  tine  and  space  scales.  The  mathematical  treatment  of  thes° 
waves  consists  in  setting  up  an  appropriate  model  for  them,  and  then  analyzina 
it.  A  natural  type  of  model  in  the  present  situation  involves  a  system  of 
partial  differential  eouations  of  reaction-diffusion  type;  small  parameters 
may  be  nut  into  the  system  to  effectuate  the  multiple  scales.  This  paper 
explains  the  basic  steps  in  modeling  chemical  waves  this  way,  and  applies  them 
to  the  case  of  taroet  patterns  (concentric  circular  waves  emanating  from  a 
point).  The  models  constructed  here  reflect  the  known  crualitative  kinetics  of 
the  px  reaction.  The  technigues  are  expected  to  be  of  value  in  reaction- 


Tbe  ">snonB  j  hi  1  i  t  v  for  t^e  wordina  and  views  exnressed  in  this  desrrintive 
summary  lies  with  Mp'’,  *nd  not  with  the  authors  of  this  report. 
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PROPAGATING  WAVES  AND  TARGET  PATTERNS  IN  CHEMICAL  SYSTEMS 


Paul  C.  Fife  and  John  Tyson 

1,  Introduction. 

The  discovery  of  propagatinq  waves  of  various  types  in  chemical  reaqents  has  provoke? 
a  great  deal  of  research,  during  the  last  ten  years,  into  the  phenomenology  and  the 
underlying  mechanisms  for  such  wavelike  activity.  The  research  has  been  performed  by 
natural  scientists  and  mathematicians  alike.  Most  of  it  has  been  experimental ,  but  much 
computer  simulation  and  mathematical  analysis  has  also  been  done.  Chemical  wave  activitv 
is  believed  to  be  prevalent  in  biological  organisms,  but  the  most  readily  accessible 
reagent  for  laboratory  study  is  that  discovered  by  Belousov  and  Zabotinskii  (the  Z- 
reagent).  This  mixture  has  oscillatory  or  excitable  kinetics,  depending  on  the 
concentrations  of  the  various  chemicals  in  the  solution.  Roth  of  these  regimes  have  at 
least  two  natural  time  scales:  During  one  period  of  an  oscillation  or  during  one  excitel 
"excursion",  most  of  the  variation  in  the  concentration  of  the  reactants  occurs  within  a 
brief  interval  of  time.  The  time  scale  associated  with  this  brief  spurt  of  activitv  is 
much  shorter  than  that  associated  with  the  slow  variation  which  occurs  before  and  after. 
This  is  well  known  from  experiment  and  computation,  and  is  evident  from  scalino  analvses  n* 
model  kinetic  equations  performed  in  [1]  and  elsewhere.  Spatial  structures  are  al«o 
prevalent  in  unstirred  layers  of  this  reagent  ([231;  [2],  (221,  [241  ,  and  references 
therein).  Target  patterns  (expanding  concentric  circular  waves)  are  amona  the  "*o«t 
prevalent  of  these  structures.  Here  again,  disparate  space  and  ti^e  scales  are  evide**- 
from  computer  simulation  of  prooagating  waves  f 3 ! . 


*This  is  a  revised  and  evpanded  version  of  a  paper  presented  ^jfn  at-  i 
conference  with  proceedings  published  as  Dynamics  of  Synergetic  Rvsfems,  . 
Haken  ,  ed ,  Rrringer-Verlag  ( 1°R0  ) . 


Sponsored  by  the  United  States  Army  under  Contract  Nos.  D  AAG2°  -7  ^  -c  -n  ^  a"  * 

PAAC,2<l-O0-C-nn41  and  National  Science  Roundation  ^ran*  No.  i, 


It  i  <5  natural,  therefor'-* ,  t.o  n<~o  m*ii  hi  r!  e  ofi’ie.j  ter  a.., , »n  r-  v.*l'’'~  ♦-/,  >-r.anr.( 

wavf>  and  pattern  phenomena  to  rn^t  l'ooa  t  ica1  analysis.  "'■•  •  *v  'r  pr^pnoT*  i  or  •lo',:)nnn  j  *j 

excitable  n»?;’ia  such  as  '•> iol o-rical  w>bnr,.o  Kr,  in  in--  :  rrof  to.1  rr  ,fTi  tht*7  uso  1  *.-■>’  .  T ^ 
general  react  ion— -1i?:  fir?  ion  set'-inqs,  those  -••*'*  tr>  ^  vore  r  ur'-qoO  i  r  ~o-.'  ■•  <  ■>  r  1  * 
ff>  ,  IP  *  Vevi_*r  the!  ,  tl’pi  •*  full  inol  ir*at  i  ops  ^a^e  ve*-  to  >  e  .{r.*  r.''i.  T  -. 

particular,  their  uj;p  in  analv/inu  tar  net  patterns  Vws  teen  ne-P  **o*-e 4 .  ",v*e  »■  irroi"  of  *-l>e 

pre.con*  paper  is  to  explain  sone  basic  principle--,  in  t  **•*  ana'!  vr:  i  s  n*  w.ivft  **•/•■>•.*«=  v*/  n-altno 
t.echnirru#»s ,  and  to  di  sciion  the  application  nt  t‘-,  ^e  t'r  i  noj  •»’.  «>•?  to  M-.-.  tin’-  or  ~»o'lel  i  no 
target  patterns.  * 

We  begin  with  a  brief  account,  in  ‘toot-ion  .7,  of  tue  f  un  i  am*  *«{-  .->  1  t^oorv  f  vr"' 
for  scalar  nonlinear  diffusion  equations,  as  ^bese  at"  f*-»  cump -north  par*  ••  -r  •■  '*«  •  •*'■ 

complex  wave  phenomena  to  b*»  examined  in  la*  or  sections.  ’  seal  ire  r»<»t*'od  *■•  *»»ee  *-'  r 

studv  of  mope  rornlcx  mrotjann*-  2  nrt  f  ronf.  s  *-q  t.'nt-  of  scalar  fronts  »  r  e  1 -»V  or in  fo'-t  i « .♦ 


3.  In  1°74,  Winfreo-  [IP  siuriest-ed  tbit-  a*  least  two  hroa 


i  O  ->  ?  V  V 


exist:  phase  anl  vM^er-  (nee  a1;;*  "}[  ).  V.’e  ran  see  qir>  j  di  s  i  on  v  •♦"  ei,*-,vi., 

in  the  context  of  wave  front'-  i:w  i  at*'  1  wi*-«  relaxation  o^c  i  1  1  i*orv ,  for  i  .»%*,•  i  t  d ]•- ) 
Vine  tics?  thic  -listinoti  -,n  i expl  a  i  mv!  too  dinrooned  at  io  . 

He sides  treating  the  i  at  ••:.:-**  mi  Properties  of  chemical  st  rur*-  ire'- ,  one  -or-  vi  r> 

inquire  how  they  mi  qb*-  arise  in  »  rea.ten*  in  rh.  •  first  pl.io*.  '”vr>  ‘-a  oh  •  e,  »-  v. ;  ip*' 

deta  i  led  in  rection  5 .  p*v«-  inn  ♦•aVes  un  *he  problem  of  n.  »•!•■•  line  *-arri»  r  1  “  i- 


f ont-i  ‘-hat-  ♦‘he  ♦,V'onrv 


. js e r . 1 1  tern.  The:"  ; 


at  lens*-  pfjme  >  f  r  i  *u»»*r  »’■ 


•  ?e.»  i  s  -nn'-Mif  i:a  1  l  * 


mi :  v  ••  i  ■  s  *  f  - 


i  ’  »  .  .*•«  f  t  •  '  irine,  •  ■  , 


Although  patterns  in  the  2-reagent  are  the  motivation  for  this  work,  specifics  of  the 
7-reagent  chemistry  are  not  touched  upon  here.  We  have  made  a  detailed  study  in  (27)  of 
the  patterned  solutions  of  a  system  of  equations  which  realistically  models  the  chemistry 
of  that  reagent. 

Other  writers  ([12,131,  and  especially  Kopell  and  Howard  fl4))  have  studied  target 
patterns  within  the  context  of  X-w  systems.  An  approach  via  Pade  approximants  is  in 
[25].  These  aporoaches  are  entirely  different  from  that  discussed  here.  The  reader 
interested  in  the  Z-reagent  can  profit  from  the  books  by  Zahotinskii  [15)  and  by  Tyson 
[16).  Some  of  the  material  in  this  paper  was  presented  from  a  different  point  of  view,  and 
in  more  mathematical  detail,  in  [17].  Some  of  the  results  announced  here  represent  ioint. 
work  with  R.  Smock.  We  wish  to  thank  M.  Marek  for  bringing  to  our  attention  papers  [4)  and 
[7). 


-1  - 


Scalar  fronts 


Here  we  review  some  basic  facts  about  wave  front  solutions  u  -  bfx-ct)  or  soa!  ir 
nonlinear  diffusion  equations 


u  -  u  +  f  ( u )  n ) 

t  xx 

where  f  has  two  zeros:  f  (n  )  =  f f tt ^ )  «*  0.  Tinder  fairly  general  circumstances ,  tber-* 
exist  fronts  satisfying 


u(-»)  =  t;  .  u(«)  =  u7. 


(? ) 


Two 


l 


important  cases  arise  in  applications: 
(a)  f*(U  )  <  0,  f'(U2)  <  ,  and  f 

a  nd  *i  ? 


has  only  one  interned iate  zero  between 


(b)  f(u)  t  0  for  u  in  the  interval  between  and  • 

In  case  (a),  there  exists  a  unique  velocity  c  and  a  profile  U(z),  unique  up  to 
shifts  in  z  /  satisfying  (2),  such  that  U(x-ct)  satisfies  (1)  [18].  This  front  is  verv 

stable;  if  it  is  perturbed  by  any  hounded  function  whose  bound  does  not  surpass  a  certain 
hnown  constant,  then  the  resulting  solution  of  (1)  evolves,  as  t  w  ,  bach  to  t1-**  sa«e 
front  { possibly  shifted  by  a  certain  amount)  (19], 

Tn  case  (b),  on  the  other  hand,  there  exists  a  whole  ranae  of  possible  front 
velocities  [20].  For  example,  suppose  f  >  0  and  >  T’0 .  Then  there  is  a  positive 

minima]  speed  c*  such  that  Cor  any  value  c,  c*  *  c  <  =°,  there  exists  a  unicine  :  *>:c-  :  ** 
for  shifts)  wave  front  solution  tJ(z),  satisfying  (2).  These  fronts  are  stable  to  -'-’nil 
perturbations  which  are  zero  except  in  a  finite  interval,  hut  they  =»re  cerf^ir.l”  not  ’ 

to  the  same  extent  as  those  of  type  (a'. 

In  anv  case,  the  front  moves  in  such  a  direction  rVy\t  for  eacv  x  , 
uf  x,  t )  =  'Ux-’t'  approaches  the  "io^inant  vr  1  or  •;  ••  a  •*  *“  - 


«!erin»»u  to  be  the  constant 
“■coral  ‘Muals  zero. 


IT 

:  {S) «'  s  <  ,  an-'  T‘.%  if  t  ii  i  s  j  at  *‘:i  ", 

1 

aep  the  f  roe*  is  c* a**  i on ar*’ . 


ir. 


then  c 


bAMU. 


3.  Decoupling  and  free  boundary  problems. 

In  typical  singular  perturbation  problems,  a  complex  system  may  be  reduced  to  several 

simpler  ones  by  rescaling  and  exploiting  the  smallness  of  some  parameter.  The  simpler 

problems  may  govern  the  solution  in  different  parts  of  its  domain  of  definition;  thus  there 

may  be  boundary  layers  versus  regions  of  relatively  slow  variation. 

Analogous  situations  arise  in  reaction-diffusion  problems  [9,10].  To  illustrate  this, 

we  consider  a  system  of  n  =  n^  +  n^  reacting  and  diffusing  components,  n^  of  them 

"fast"  and  the  others  not  fast.  Let  u  «  (u„ ,...,u  )  be  the  vector  of  fast  components, 

1  n1 

and  v  the  others.  The  system  of  RD  equations  in  one  space  variable  is  of  the  form 


aD.u  +  kf(u,v). 

(3a) 

1  XX 

=  D,v  +  g(u,v). 

2  xx 

(3b) 

Here  are  (diffusion)  matrices;  k  >>  1  is  a  parameter  expressing  the  fact  that 

reactions  affecting  the  concentration  u  are  fast;  and  the  parameter  a  <<  k  is  inserted 
to  account  for  the  possibility  that  the  diffusion  rate  of  u  may  be  small  or  large 
(D^  =0(1))  .  The  lowest  order  approximation,  in  regions  where  ut  and  ciu  ^  are  not 

large,  is  obtained  by  setting  the  coefficient  of  k  equal  to  zero: 

f  ( u ,  v )  =  0.  (4) 

We  assume  that  this  equation  can  be  solved  for  u  in  a  nonunicue  manner:  there  are  at 

n2  n, 

least  two  functions  h+ ,h  :  R  ♦  R  ,  such  that  (4)  holds  when 

u  =  h±(v) .  (5  ) 

Tn  other  words,  f(h+(v),v)  5  f(h  (v),v)  =  0. 

We  imagine  that  the  x-t  plane  is  partitioned  into  two  parts  in  which  (5)  holds 

(approximately)  with  the  corresponding  sign.  Under  certain  circumstances,  such  a 
partitioning  is  possible,  with  sharp  wave  fronts  forminq  the  boundaries  between  the  two 
domains. 

To  investigate  this  further,  we  scale  x  and  t  differently  in  the  layer  ''etwee-, 
w  and  i!  .  The  scaling  will  be  chosen  so  as  to  eliminate  the  parameters  in  (3a). 


MM 


In  more  general  situations,  the  boundary  between  and  il  could  consist  of 
several  wave  fronts.  And,  of  course,  the  extension  of  this  reasoning  to  higher  space 
dimensions  is  clear. 
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4  Phase  and  trigger  fronts. 

We  refer  to  the  setting  in  Section  3 
higher  dimensions  remain  to  be  explored, 
additional  complication  we  impose  is  that 
rather  their  graphs  lie  on  a  nullcline  of 


with  n1  =  n^=  1;  corresponding  phenomena  is 
So  now  a  and  v  are  scalar  functions.  The 
h^  are  not  defined  for  all  values  of  v  ; 
f  as  shown: 


Figure  1 


For  values  of  v  in  the  interval  v_  <  v  <  v  ,  f(u,v)  has  the  features  of  the  function  i 

(1)  in  case  (a),  provided  f  (h,(v),v)  *  0. 

u  ± 

'■!e  assume,  mainly  for  simplicity,  that  =  0  in  (3b),  and  that 

n(h_(vl,v)  _  r^(v)  <  0,  g(h  (v),v)  =  G?(v)  >  0.  If  there  is  a  sinole  front  wit'" 
trajectory  y(t)  and  v  varying  continuously  across  it,  then 


3v 

3t 


c,_(v)  <  n  ,  x  <  v(t)  , 
G+  ( v )  >0  ,  x  >  v  f  t )  . 


If  v  <  v  <  V  ,  y(t)  is  ooverned  by  (7),  because  the  function  O'"' 


in  t'  t- 


-n-n.i-  consideration  of  the  bistable  case  ((a),  Fortier  g).  Further  -ore ,  ■  c 


continuous  at  fho  front  (as  we  assume' 


if  must  remain  so 


For 


front  passes  a  fixed  value  of  x  ,  v  would  change  discontinuously  in  time,  meaning  than 


on 


v  ).  In  this  sense,  the  motion  of  phase  fronts  is  determined  by  the  ini 


branches,  say  h  .  After  this  happens,  a  slower  process  takes  place,  in  which  v  evolves 
according  to  (3b)  with  u  replaced  by  h_(v).  Again  for  simplicity,  assume 
D,  =  0,  so  v  =  G_(v)  <  0.  Eventually,  for  some  x  =  xQ,  v  will  attain  a  minimal 
value  of  v  •  Further  decrease  of  v  causes  the  image  to  leave  the  branch  u  =  h  (v) 
for  x  in  a  neighborhood  of  x^.  Then  by  the  process  described  in  (i)  above,  that  part 
of  the  image  curve  is  rapidly  attracted  to  the  other  stable  branch  u  =  h+(v).  This 
localized  attraction  to  h+  causes  a  pair  of  fronts,  facing  oppositely,  to  be  formed  near 
x  =  xQ  .  For  each  such  front,  h+  will  be  the  dominant  state,  so  the  fronts  will  move 
apart,  increasing  the  interval  on  which  u  ~  h+(v). 


These  are  a  series  of  concentric  circular  chemical  waves,  expanding  outward,  new  ores 
regularly  being  generated  at  the  center  (usually  called  a  "leading  center"  [21]).  Such 
patterns  have  been  observed  in  various  forms  of  the  Z-reagent  ([23);  [2],  [22],  and 
references  therein).  Some  of  the  patterns  observed  are  associated  with  externally  imposed 
heterogeneities  at  the  center.  We  shall  indicate  in  (i)  below  how  such  targets  may  be 
modelled.  Self-sustaining  target  patterns,  not  dependent  upon  external  stimuli,  may  also 
be  modelled  by  the  techniques  discussed  above;  see  (ii)  below. 

All  of  the  models  we  describe  involve  the  same  two  basic  phenomena:  (a)  spontaneous 
generation  of  wave  fronts  at  the  leading  centers  as  described  in  Section  5  (ii),  and  (b) 
their  subsequent  motion,  according  to  the  rules  brought  out  in  Section  4.  The  generation 
process  in  5  (ii)  was  for  pairs  of  diverging  fronts  moving  in  one  space  dimension;  its  two- 
dimensional  analog  is  the  spontaneous  appearance  of  a  small  circular  front  which  spreads 
outward.  Since  the  fronts  are  very  narrow  in  our  analysis,  they  appear  locally  as  plane 
waves.  Therefore  it  suffices  to  treat  the  problem  in  a  one-dimensional  framework,  which  we 
shall  do.  Alternatively,  the  variable  x  could  be  interpreted  as  distance  to  the  origin 
in  a  configuration  with  radial  symmetry. 

(i)  Imposed  heterogeneities.  We  suppose  that  near  the  origin  there  is  a  substance 
with  prescribed  density  distribution  w(x);  alternately,  w  could  represent  an  imposed 
temperature  distribution.  We  also  suppose  that  w  influences  the  reaction  process,  so 
that  f  and  y  in  (3)  are  functions  of  u,v,  and  w.  For  each  value  of  w  ,  the 
nullcurves  f  *  0  and  g  =  0  have  the  shape  shown  in  Figure  1  ;  but  their  relative 
positions  may  vary  with  w  .  In  particular,  h+,  v  ,  and  v  may  depend  on  w  . 

At  x  =  0,  fronts  involving  an  abrupt  increase  in  u  (upjump  fronts)  form,  no  -  .r  : .  r  • 
to  the  description  in  Section  5,  when  v  has  a  local  minimum  at  the  origin  wMc'* 
to  v  .  Similarly,  downjump  fronts  form  when  u  =  h+(v)  and  v  has  a  maximum  whic'- 
increases  to  v  .  The  result  in  that  when  x  is  fixed  at  0  ,  the  trajectory 
( u(  n  ,*■ ) ,  v(  h ,  t ) )  is  that  of  a  relaxation  oscillator  with  the  kinetics  o'  riour.-  1.  s,  • 
nullcurve  g  -  n  must  be  placed  as  shown,  to  ensure  that  the  kinetic;  ire  osci 1 )  o  a--. 


The  pattern  must  be  periodic  in  time.  The  period  T  is  set  by  the  period  of  the 
relaxation  oscillatory  motion  followed  by  the  solution  at  x  =  0,  w  =  w(0).  This  is 
presumed  known. 

The  mathematical  analysis  of  the  above  conceptual  model  consists  in  determining  the 
trajectories  of  the  expanding  fronts,  and  the  function  v(x,t),  so  that  all  the  above 
evolutionary  laws  and  constraints  are  fulfilled.  It  is  convenient  to  express  the  fronts' 

j  + 

motion  by  the  functions  T  (x).  Here  T  (x)  is  the  time  at  which  some  specific  up jump 
front  reaches  position  x,  and  t  (x)  +  is  the  time  for  the  next  succeeding  downjump.  The 
next  upjump  is  then  at  time  T  +  T  .  The  equations  to  be  satisfied  are: 

t±(x)  «  t  c  1(v(x,T±(x) ) ,w(x) ) ,  (for  trigger  fronts), 


d_ 

dx 


T±(x) 


given  by  (9)  (for  phase  fronts). 


3v 

3t 


f 

v. 


G+( v,w(x> ) ,T+( x)  <  t  <  T  (x)  , 

G  (v,w< x) ) ,  T  (x)  <  t  <  T  (x)  +  T 


+  + 

(it  suffices  to  determine  v  in  the  interval  T  <  t  <  x  +  T  ).  The  periodicity 

+  + 

constraint  on  v  is  that  v(x,x  ( v ) )  =  v(x,T  (x)  +  T).  In  addition,  we  must  require 
that  v(0,t)  be  the  values  of  v  corresponding  to  the  relaxation  oscillator  at  the 
center,  and  that  as  x  ■*  “>  ,  the  x*(x)  approach  linear  functions  (corresponding  to  a 
plane  wave  train). 

If  we  insist  that  all  the  fronts  be  trigger,  then  the  ineauality 

v(w(x) )  <  v  <  v(w(x) ) 

is  an  additional  constraint,  and  it  is  a  nontrivial  matter  to  determine  whether  there  exist 
functions  v  ,  satisfying  all  the  above.  However,  the  problem  becomes  somewhat  nasi—- 
one  allows  every  other  front  (say,  all  the  down  jump  fronts'  to  be  of  -Viaso  tvre,  ‘•'•e 
remaining  triqqer.  This  appears  to  be  the  situation  common)  ••  observe  !  in  *ho  "-rein,?n‘ 
ta  rgets . 
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Tvson  1  ha*  *o  '>■  1 1  n  *■  *  scaled  version  of  *  'i*r  of  react  ion-di  rf  usi  on  equations  (the 

"Oregonattr" '  realistically  node  lino  the  Belousov-habot i nsk i  1  reaction.  We  [27]  have 

invest  i  natod  turrets  with  al-^r-iatine  phase  and  trigger  fronts  usina  this  system. 

The  third  conceptual  possibility  is  when  alL  fronts  are  of  phase  type.  Then  the 
conf  i  omr  at  ion  of  fie  rr-v>t\s,  as  **hoy  *ev'»lor>  in  the  reaoont,  will  depend  completely  on  the 

initial  concent rat io-s ,  and  will  not  in  general  be  circles.  Such  patterns,  if  they  do 

exist,  would  therefore  not  account  eor  the  observed  circular  fronts. 

fii)  Fe 1 f-susta i n i na  tarn*=*r.  patterns  may  be  modelled  by  retain  inn  the  function  w(x), 
but  supposing  to  obev  a  thi^d  react ion-di f fusion  equation  coupled  to  the  first  two. 

Thus  the  distribution  of  w  i<  obtained  as  part  of  the  model,  not  imposed  by  external 
conditions.  TMs  type  of  mo  lei  is  explored  in  some  detail  in  [17],  assuming  all  fronts  are 
trigger.  Again,  the  dynamics  of  the  wave  front  and  the  functions  v  and  w  involve 
complicated  mathematics,  but  s.  me  reasonable  simplifications  are  possible. 

A  conceptually  similar  model  for  self-susta ininq  patterns  involving  three  reacting 
components  had  previously  been  proposed  by  Zaihin  and  Kawczvnski  [28];  but  its  analysis 
was  left  incomplete  and  question  of  phase  and  trigger  fronts  was  not  addressed. 


7.  Targets  in  excisable  media 


With  nullcurves  as  depicted  in  Fiqure  1, 

ufc  =  kf(u,v) 

have  stable  relaxation  oscillatory  solutions, 
the  positions  in  Figures  3  and  4,  however,  the 


vf 

l 


Figure  3 


the  kinetic  eauations 
v  =  g(  u ,  v) 

When  the  g  nullcurve  is  shifted  to  one  of 
kinetics  become  excitable  or  bistable. 

v  i 


Figure  4 


For  example  in  Fiqure  3,  there  is  one  stable  rest  state  (the  uniaue  intersection  point); 
but  when  this  state  is  perturbed  downward  a  small  amount,  the  solution  makes  a  laroe 
excursion  (dotted  line)  before  returning  to  the  rest  state. 

In  modeling  target  patterns  in  Section  6,  we  needed  the  kinetics  to  be  oscillatory  at 
the  oriqin;  but  away  from  the  origin,  the  conf igurations  in  Fiqures  3  and  4  are  nor 
excluded.  As  x  *  00  the  tarqet  develops  into  a  reqular  wave  train,  so  of  course  the 
kinetics  must  support  such  a  train  there.  Excitable  and  bistable  (as  well  as  oscillatory) 
kinetics  do  support  wave  trains.  For  example,  the  phase  plane  imaqe  (orbit)  of  a  train  it 
a  fixed  value  of  x  is  shown  by  the  dotted  loon  in  Figure  4.  Therefore  ther^  is  no 
contradiction  involved  in  having  a  periodic  target,  pattern  emerge  in  an  excitable  medium. 
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This  is  apparently  what  often  happens  with  the  Z-reagent,  which  can  exist  in  an 
excitable/  as  well  as  oscillatory,  regime.  The  initial  formation  of  a  target  pattern  in 
such  a  medium  involves  a  wave  train  entering  a  guiescant  region  at  a  stable  rest  state.  It 
may  be  difficult  to  visualize  this  process  of  excitation  into  a  periodic  state,  but  several 
mechanisms  are  available  for  accomplishing  it  [24;  R.  Smock,  in  preparation] . 


-16- 


References 

1.  J.  Tyson  1979,  Oscillations,  bistability,  and  echo  waves  in  models  of  the  Belousov- 
Sfhabotinskil  reaction.  Ann.  N.  Y.  Acad.  Sci.  36,  279-295. 

2.  A.  T.  Winfree  1978,  Stably  rotating  patterns  of  reaction  and  diffusion, 

Theor.  Chem.,  Vol.  4,  Academic  Press,  New  York,  pp.  1-51. 

3.  D.  J.  Reusser  and  R.  J.  Field  1979,  The  transition  from  phase  waves  to 
trigger  waves  in  a  model  of  the  Zhabotinskii  reaction.  J.  Amer.  Chem. 

Soc.  101,  1063-1071. 

4.  L.  A.  Ostrovskii  and  V.  G.  Yakhno  1975,  The  formation  of  pulses  in  an 
excitable  medium.  Biofizika  20,  489-493. 

5.  R.  Casten,  H.  Cohen,  and  P.  Lagerstrom  1975,  Perturbation  analysis  of  an 
approximation  to  Hodgk in-Huxley  theory.  Quart.  Appl.  Math.  32,  365-402. 

6.  J.  P.  Keener  1979,  Waves  in  excitable  media,  preprint. 

7.  V.  G.  Yakhno  1975,  On  a  model  for  leading  centers,  Biofizika  20,  669-674. 

See  also  G.  M.  Zhislin,  V.  G.  Yakhno,  and  Ju.  X.  Goltsova  1976,  Biofizika 
21,  692-697;  Ju.  K.  Goltsova,  G.  M.  Zhislin,  and  V.  G.  Yakhno  1976, 

Biofizika  21,  893-897;  V.  G.  Yakhno  1977,  Biofizika  22,  876-881. 

8.  G.  Carpenter  1978,  Bursting  phenomena  in  excitable  membranes,  SIAM.  J. 

Appl.  Math.,  to  appear. 

9.  P.  C.  Fife  1976,  Pattern  formation  in  reacting  and  diffusing  systems,  J. 

Chem.  Phys.  64,  854-864. 

10.  p.  c.  Fife  1976,  Singular  perturbation  and  wave  front  techniques  in 
reaction-diffusion  problems,  in:  SIAM-AMS  Proceedings,  Symposium  on 
Asymptotic  Methods  and  Singular  Perturbations,  New  York,  23-49. 

11.  A.  T.  Winfree  1974,  Wavelike  activity  in  biological  and  chemical  media, 
in:  Lecture  Notes  in  Biomathematics  (Ed.  P.  van  den  nriessche), 

Springer-Verlag,  Berlin. 

12.  J.  M.  Greenberg  1978,  Axisymmetric  time-periodic  solutions  of  reaction-diffusion 
equations,  SIAM  J.  Appl.  Math.  34,  391-397. 


-17- 


13.  P.  Ortoleva  dP.  1  J.  w :>ss  1974,  On  a  variety  of  wave  phenomena  in  chemical 
and  biochemical  rmrilUtions,  J.  Chem.  Phys.  60,  5000-5107, 

14*  N’.  Kopoll  and  L.  \*.  Howard  11'79.  Target  patterns  anii  horseshoes  from  a 

pert.:rbei  central  force  problem:  some  temporally  periodic  solutions  to 
reao?  o'**'  diffusion  equations.  Preprint. 

If.  A.  ".  Oabotinskn  1974,  Concentration  Oscillations  (Russian),  Nauka,  Moscow. 

[•- .  J.  J.  ’TV so»'.  i°7f»,  The  Belnusov-Zhabot inski i  Reaction,  Lecture  Notes  in 
Bicvnether.at  ies  No.  10,  Springer,  New  York. 

17,  ?,  c.  Fife  1^70,  Wave  fronts  and  target  patterns,  in:  Applications  of 

Nonlinear  Analysis  in  the  Physical  Sciences,  Pitman  Publishing,  London, 
to  appear. 

IS.  Ya.  I.  Kanel*  1962,  On  the  stabilization  of  solutions  of  the  Cauchy 
problem  for  the  equations  arising  in  the  theory  of  combustion.  Mat. 

Sbornik  59,  245-288. 

19.  P.  C ,  Fife  and  J.  B.  McLeod  1977,  The  approach  of  solutions  of  nonlinear 
diffusion  equations  to  traveling  front  solutions.  Arch.  Rational  Mech. 

Anal.  65,  335-361.  Also:  Bull.  Amer.  Math.  Soc.  81,  1075-1078  (1975). 

20.  A.  N.  Kolmogorov,  I.  G.  Petrov' skix,  and  N.  S.  Piskunov  1937,  A  study  of 
the  equation  of  diffusion  with  increase  in  the  quantity  of  matter,  and  its 
application  to  a  biological  problem,  Bjul.  Moskovskovo  Cos.  Univ.  17,  1-72 

21.  A.  M.  Zhabotinsky  and  A.  N.  Zaikin  1973,  Autowave  processes  in  a  distributed 
chemical  system,  J.  Theor.  Biol.  40,  4S-61. 

22.  M.  Marek  and  J,  Juda  197",  Controlled  generation  of  react ion-dif f usion  waves,  Sci. 
papers  of  Prague  Inst.  Che*".  Technol.  S °r.  Y.,  to  appear. 

23.  n.  7a  iV  \  n  and  A.  *•*..  Zhabot insky  1^70,  Concentration  wave  propagation  in  two- 
j  j  rr»ecs  iocal  \  i  mid-phase  self-oscillating  system,  Nature  225,  535-537. 


-  1  '•»  - 


24.  M.  -L.  Smoes  1990,  Chemical  waves  in  the  oscillatory  ^Jhabotinskii  system;  a  transition 

from  temporal  to  spatio-temporal  organization,  pp.  80-96  in  Dynamics  of  Synergetic 
Systems,  H.  Haken,  ed.,  Springer  Verlag,  Berlin. 

25.  P.  Ortoleva  1978,  Dynamic  Pade  approximants  in  the  theory  of  periodic  and  chaotic 
chemical  center  waves,  J.  Chem.  Phys.  69,  300-307. 

26.  P.  Ortoleva  and  J.  Ross  1975,  Theory  of  propagation  of  discontinuities  in  kinetic 
systems  with  multiple  time  scales:  fronts,  front  multiplicity,  and  pulses,  J.  Chem. 
Phys.  63,  3398-3408. 

27.  J.  Tyson  and  P.  C.  Fife  1980,  A  realistic  model  for  target  patterns  in  the  Belousov- 
Xhabotinskii  reaction,  to  appear. 

28.  A.  N.  2aikin  and  A.  L.  Kawczynski  1977,  Spatial  effects  in  active  chemical  systems. 

I.  Model  of  leading  center,  J.  Non-Equilib.  Thermodyn.  2,  39-48. 


PCF/JT/db 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 

2074 


4.  TITLE  (end  Subtitle) 


PROPAGATING  WAVES  AND  TARGET  PATTERNS  IN  CHEMICAL 
SYSTENS 


7.  AUTHORC«J 

Paul  C.  Fife  and  John  Tyson 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


t».  CONTROLLING' OFFICE  NAME  AND  ADDRESS 

See  18  below 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  A  PERIOD  COVEREO 

Summary  Report  -  no  specific 
reporting  period 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERS e) 

DAAG29-80-C-0041 
DAAG29-7  5-C-0024 
MCS79-04443 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  8  WORK  UNIT  NUMBERS 

2  -  Physical  Mathematics 


IZ.  REPORT  DATE 

May  1980 


13-  NUMBER  OF  PAGES 


I  X  L' 


4.  MONITORING  VGENCY  NAME  4  ADDRESS, 'if  dilierent  from  Controlling  Olfice)  IS.  SECURITY  CLASS,  (of  (hi*  fport) 

UNCLASSIFIED 


16.  OISTP  BUTION  STATEMENT  (ot  thle  Report) 

Approved  for  public  release;  distribution  unlimited. 


17-  DISTRIBUTION  STATEMENT  (ot  the  ebetract  entered  In  Block  20,  II  dlllerent  trom  Report) 


18.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  office  National  Science  Foundation 

P.O.  Box  12211  Washington,  D.C.  20550 

Research  Triangle  Park, 

North  Carolina  27709 


19.  KEY  WOROS  (Continue  on  r«v«ra«  aid*  if  naceaamry  and  identify  by  block  number) 

Reaction,  Diffusion,  Chemical  waves,  Wave  front,  Pattern  formation,  Singular 
perturbations 


JO  ABSTRACT  (Continue  on  reverse  aide  II  necesaary  end  Identity  by  block  number) 

A  class  of  models  for  target  patterns  (concentric  circular  waves  emanating 
from  a  point  called  the  leading  center)  is  constructed  in  the  context  of  singu¬ 
larly  perturbed  reaction-diffusion  systems  of  partial  differential  equations. 
First,  the  theory  of  wave  fronts  is  detailed  for  scalar  equations  and  systems  of 
equations.  A  scaling  method  reduces  complex  waves  to  the  consideration  of  a 
group  of  simple  wave  phenomena.  It  is  shown  that  expanding  wave  fronts  can  be 
generated  spontaneously  at  a  point.  This  process,  together  with  the  laws  of 
their  subsequent  motion,  reduces  the  problem  to  an  ordinary  differential  initial 


EDITION  OF  I  NOV  88  IS  OBSOLETE 


UNCLASSIFIED  (continued) 

security  CLASSIFICATION  OF  TH  S  PAGE  ('Wian  Date  Bntered) 


20  Abstract  (continued) 

value  problem,  whose  solution  is  required  to  have  certain  properties.  A  discussion 
is  given  of  the  connection  between  these  results  and  experimental  observations. 


